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The  object  of  this  report  Is  to  derive  the  general  hydro- 
dynamical  theory  needed  for  mathematical  treatment  of  problems 
involving  the  Interaction  between  bodies  immersed  in  water  and 
surface  gravity  waves  in  the  water.   In  particular,  the  problem 
of  the  effectiveness  of  floating  breakwaters  in  the  form  of 
rigid  bodies  or  of  beams  with  bending  flexibility  will  be 
treated;  In  these  cases  the  principal  object  is  to  determine 
the  amount  of  reflection  of  an  Incoming  progressing  wave  that 
can  be  achieved  by  such  structures  in  terms  of  parameters  which 
characterize  the  wave  motion  and  the  physical  nature  of  the 
floating  breakwater. 

The  theory  presented  here  is  an  approximate  theory  derived 
from  the  exact  hydrodynamical  theory  for  gravity  waves  of  small 
surface  amplitude  by  making  the  simplifying  assumption  that  the 
wave  length  X  in  the  water  is  sufficiently  large  compared  with 
the  depth  h  of  the  water.   If  the  ratio  X/h  is  of  the  order 
eight  to  ten  the  results  should  be  reasonably  accurate,  and 
such  a  condition  would  often  be  realized  in  cases  having  prac- 
tical interest.   It  is  entirely  possible  to  formulate  a  mathe- 
matical theory  without  making  such  a  restrictive  assumption 
regarding  the  ratio  of  wave  length  to  depth  —  for  floating 
rigid  bodies  having  simple  harmonic  motions,  for  example,  the 
theory  has  been  worked  out  by  P.  John  --  but  the  difficulties 
involved  in  obtaining  solutions  in  cases  of  Interest  In  practice 
are  so  formidable  that  no  progress  In  this  direction  has  been 
achieved. 

The  material  is  presented  here  in  three  parts.   The  first 
section  (by  J.  J.  Stoker)  consists  of  a  derivation  of  a  theory 
of  floating  bodies  In  shallow  water  which  is  carried  out  in 
such  a  way  as  to  permit  floating  bodies  of  widely  different 
physical  characteristics:  rigid  bodies,  inert  distributions  of 
mass,  membranes  under  tension,  or  beams  with  bending  stiffness, 
for  example.   In  Section  2  (by  B.  Fleishman),  the  details  of 
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the  theory  are  worked  out  for  floating  rigid  bodies  and  numerical 
results  are  given  for  a  floating  board  having  one  degree  of 
freedom  (a  purely  vertical  motion),  and  two  degrees  of  freedom 
(angular  and  vertical  motion).   Reflection  and  transmission 
coefficients  are  calculated  as  well  as  the  motion  of  the  board 
and  the  dynamic  pressure  on  it.   The  results  are  shown  in  the 
form  of  curves.   In  Section  3  (by  L.  Weliczkor)  the  relevant 
theory  is  derived  and  applied  to  floating  beams  of  infinite, 
semi-infinite,  and  finite  length.   (It  might  be  mentioned  that 
this  was  the  most  difficult  part  of  the  investigation.) 

The  case  of  a  floating  beam  was  suggested  to  us  by   J.  H. 
Carr  of  the  Hydraulics  Structures  Laboratory  at  the  California 
Institute  of  Technolofy  as  one  having  practical  possibilities; 
at  his  suggestion  calculations  in  specific  numerical  cases  have 
been  carried  out  in  order  to  determine  the  effectiveness  of 
such  a  breakwater.   The  reason  for  considering  such  a  structure 
for  a  breakwater  as  a  means  of  creating  relatively  calm  vjater 
between  it  and  the  shore  Is  the  following:  a  structure  which 
floats  on  the  surface  without  sinking  far  into  the  water  need 
not  be  subjected  to  large  horizontal  forces  and  hence  would  not 
necessarily  require  a  massive  anchorageo   However,  in  order  to 
be  effective  as  a  reflector  of  waves  such  a  floating  structure 
would  probably  have  to  be  built  with  a  fairly  large  dimension 
in  the  direction  of  travel  of  the  Incoming  waves.   As  a  conse- 
quence of  the  length  of  the  structure,  it  v.'ould  be  bent  like  a 
beam  under  the  action  of  the  waves  and  hence  could  not  in 
general  be  treated  with  accuracy  as  a  rigid  body  In  determining 
its  effectiveness  as  a  barrier.   This  brings  with  it  the  pos- 
sibility that  the  structure  might  be  bent  sufficiently  that  the 
stresses  set  up  would  be  a  limiting  feature  in  the  design.   The 
main  object  of  the  present  report  Is  to  present  a  means  of 
studying  all  of  these  questions  thev-)retlcally  and  to  make  cal- 
culations in  a  number  of  specific  cases. 

The  floating  beam  is  assumed  to  be  placed  at  right  angles 
to  the  direction  of  the  incoming  waves  (i.e.  its  long  axis  is 
at  right  angles  to  the  wave  crests).   The  water  is  supposed  to 
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extend  to  Infinity  In  both  of  the  directions  along  the  axis  of 
the  beam.   The  specifications  (as  suggested  by  Carr)  for  a  beam 
having  a  width  of  one  foot  (parallel  to  the  wave  crests,  that  Is) 
were  : 

Weight :   85  pounds/foot 

Moment  of  Inertia  (of  area  of  cross-section):   0.2(ft,) 

6  2 

Modulus  of  elasticity:  30x10  (pounds/in.  )   or 

U37x  lo'^(pounds/ft.)^ 

The  depth  of  the  water  Is  taken  as  [|.0  feet.   Simple  harmonic 
progressing  waves  having  periods  of  6  and  of  l5  sees,  were  to 
be  considered,  and  these  correspond  to  wave  lengths  of  250  and 
507  feet,  and  to  circular  frequencies  cr  of  l|l8  x  lo"   and 
785x10   cycles  per  second,  respectively.   The  problem  Is  to 
determine  the  reflecting  power  of  the  beam  under  these  circum- 
stances when  the  length  of  the  beam  is  varied.   In  other  words, 
we  assiime  a  wave  train  to  come  from,  say,  the  right  hand  side 
of  the  beam  and  that  it  Is  partly  transmitted  under  the  beam  to 
the  left  hand  side  and  partly  reflected  back  to  the  right  hand 
side.   The  ratio  R/B  of  the  amplitude  R  of  the  reflected  wave 
and  the  amplitude  B  of  the  incoming  wave  is  a  measure  of  the 
effectiveness  of  the  beam  as  a  breakwater.   The  effectiveness 
of  the  breakwater  can  also  be  measured  in  terms  of  the  trans- 
mitted and  reflected  energies.   If  T  is  the  amplitude  of  the 

transmitted  wave  the  energy  balance  Is  expressed  by  the  relation 

2         2  '2 

(T/B)   +  (R/B)   =  1,  so  that  (R/3)  ,  for  example,  represents 

the  fraction  of  the  Incoming  energy  which  is  reflected. 

Before  discussing  the  case  of  beams  of  finite  length  it  is 
interesting  and  worthwhile  to  consider  seml-lnflnite  beams 
first.   Since  the  calculations  are  easier  than  for  beams  of 
finite  length  It  was  found  possible  to  consider  a  larger  range 
of  values  of  the  parameters  than  was  given  above.   The  results 
are  summarized  in  the  following  tables: 
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Table  A 


X  (ft) 

^(sec) 

W( pounds) 

I  (rt)^ 

^   /pounds\ 
V    ft        / 

h   (ft) 

F./B 

507 

.lj.18 

85 

.20 

I;37xlo''' 

ko 

.li+ 

250 

.765 

85 

.20 

1^37x10''' 

ko 

.19 

196 

1.0 

65 

.20 

k31^  lo'' 

hQ 

.23 

114-7 

1.5 

85 

.20 

i|37x  lo'' 

~  ko 

.32 

98 

2.0 

85 

.20 

i|37xlO' 

ko 

.43 

In  this  table  the  beam  design  data  are  as  given  above.   At  the 
two  specified  circular  frequencies  of  O.I4.I8  and  0.765  one  sees 
that  the  floating  beam  Is  quite  Ineffective  as  a  breakwater  since 
the  reflected  wave  has  an  amplitude  of  less  than  1/5  of  the 
amplitude  of  the  Incoming  wave,  even  for  the  higher  frequency 
(and  hence  shorter  wave  length),  which  means  that  less  than  k  /o 
of  the  Incoming  energy  Is  reflected  back.   At  higher  frequencies, 
and  hence  smaller  wave  lengths,  the  breakwater  Is  more  effective, 
as  one  would  expect.   However,  the  approximate  theory  used  to 
calculate  the  reflection  coefficient  R/B  can  be  expected  to  be 
accurate  only  If  the  ratio  X/h  of  wave  length  to  depth  Is  suf- 
ficiently large,  and  even  for  the  case   X  =  25o  ft.  ( o~  =  .765) 
the  reflection  coefficient  of  value  0.19  may  be  in  error  by 
perhaps  lO"/)  or  more  since  X/h  Is  only  about  6,  and  the  errors 
for  the  shorter  wave  lengths  would  be  greater. 

Tatle  B 


X 

C" 

w 

I 

E 

h 

R/B 

507 

.418 

384 

o20 

437x10"^ 

ko 

.51 

250 

.785 

3614- 

.20 

437x10"^ 

ko 

.75 

The  only  change  as  compared  with  the  first  two  rows  of  table  A 
la  that  the  weight  per  foot  of  the  beam  has  been  Increased  by  a 
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factor  of  more  than  I|.  from  a  value  of  85  pounds/ft  to  a  value 
of  3814.  pounds/ft.   The  result  Is  a  decided  increase  In  the 
effectiveness  of  the  breakwater,  especially  at  the  shorter  wave 
length,  since  more  than  half  (i.e.  (.75)  )  of  the  incoming 
energy  would  be  reflected  back.   However,  this  beneficial  effect 
is  coupled  with  a  decided  disadvantage,  since  quadrupling  the 
weight  of  the  beam  would  cause  it  to  sink  deeper  in  the  water 
in  like  proportion  and  hence  might  make  heavy  anchorages  ne- 
cessary. 

Table  G 


\ 

cr 

w 

I 

E 

h 

R/B 

507 

.[^18 

85 

2.0 

14.37  X  10 "^ 

ko 

.26 

250 

.785 

65 

2.0 

i|37x  lo"^ 

ko 

.32 

CO 

1 

Table  C  Is  the  same  as  the  first  two  rows  of  table  A  except  that 
the  bending  stiffness  has  been  Increased  by  a  factor  of  10  by 
increasing  the  moment  of  inertia  of  the  beam  cross-section  from 
0.2  ft^  to  2.0  ft^.   Such  an  Increase  in  stiffness  results  in  a 
noticeable  Increase  in  the  effectiveness  of  the  breakwater,  but 
by  far  not  as  great  an  increase  as  is  achieved  by  multiplying 
the  weight  by  a  factor  of  four.   If  the  stiffness  were  to  be 
made  infinite  (i.e.  if  the  beam  were  made  rigid)  the  reflection 
coefficient  could  be  made  unity,  and  no  wave  motion  would  be 
transmitted.   This  is  evidently  true  for  a  semi-Infinite  beam, 
but  would  not  be  true  for  a  rigid  body  of  finite  length. 

Table  D 


X 

o~ 

A 

T 

E 

h 

R/B 

507 

.[4.18 

85 

0 

0 

ko 

.  001 

250 

.765 

85 

0 

0 

ko 

.007 
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Here,  the  difference  as  compared  with  table  A  is  that  the  beam 
stiffness  Is  taken  to  be  zero.   This  means  that  the  surface  of 
the  water  Is  assumed  to  be  covered  by  a  distribution  of  inert 
particles  weighing  85  pounds  per  foot.   As  we  observe,  there  Is 
practically  no  reflection  and  this  is  not  surprising  since  the 
mass  distribution  has  a  value  that  Is  only  slightly  greater 
than  that  of  water  Itself  (62. 1|  pounds  per  foot). 

One  might  perhaps  summarize  the  above  results  as  follows: 
A  very  long  beam  can  be  effective  as  a  floating  breakwater  if 
it  is  stiff  enough.   However,  a  reasonable  value  for  the  stiff- 
ness (the  value  0.2  given  above)  leads  to  an  ineffective  break- 
water unless  the  specific  weight  of  the  beam  is  fairly  large 
multiple  (say  8  or  10)  of  the  weight  of  water. 

In  practice  It  seems  unlikely  that  beams  long  enough  to 
be  considered  semi-infinite  would  be  practicable  as  breakwaters. 
(The  term  "long  enough"  might  be  Interpreted  to  mean  a  sufficiently 
large  multiple  of  the  wave  length,  but  since  the  wave  lengths 
are  of  the  order  of  200  feet  or  mere  the  correctness  of  this 
statement  saens  obvious.)   It  therefore  is  necessary  to  inves- 
tigate the  effectiveness  of  beams  of  finite  le?:gth.   Such  an 
investigation  requires  extremely  tedious  calculations  —  so  much 
so  that  only  a  certain  number  of  niomerical  cases  have  been  treated. 
These  are  summarized  in  the  following  tables. 


0-=  .785 


X  =  250 


X  =  507 


(ft) 

R/3 

17.5 

0 

i+9-2 

.93 

72.9 

0 

98.5 

.75 

11+5  o  9 

.10 

196.9 

0 

291.8 

.33 

U5o.i^ 

>32 

583.6 

.12 

656.3 

.13 

875.1; 

.32 

oo 

.19 
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In  these  tables  the  parameters  have  values  the  same  as  in  the 
first  two  rows  of  table  A,  except  that  now  lengths  other  than 
infinite  length  are  considered.   The  most  noticeable  feature  of 
the  results  given  in  the  tables  is  their  irregularity  and  the 
fact  that  at  certain  lengths — even  certain  rather  short  lengths-- 
the  beam  proposed  by  Carr  seems  to  be  quite  effective.   For 
example,  when  the  wave  length  is  200  ft,  a  beam  less  than  $0   ft. 
long  reflects  more  than  8o°/o  of  the  incoming  energy,   A  beam  of 
length  I4.U3  ft.  is  also  equally  effective  at  the  longer  wave 
length  of  507  ft."^ 

However,  the  maximiom  effectiveness  of  any  such  breakwater 
occurs  for  a  specific  wave  length  within  a  certain  range  of  wave 
lengths;  thus  the  reflection  of  a  given  percentage  of  the  in- 
coming wave  energy  would  involve  changing  the  length  (or  some 
other  parameter)  of  the  structure  in  accordance  with  changes  in 
the  wave  length  of  the  incoming  waves.   Also,  the  reflection 
coefficient  seems  to  be  rather  sensitive  to  changes  in  the 
parameters,  particularly  for  the  shorter  structures:   a  rela- 
tively slight  change  in  length  from  an  optimum  value  (or  a 
slight  change  in  frequency)  leads  to  a  sharp  decrease  in  the 
reflection  coefficient.   It  is  also  possible--even  probable — 
that  the  shallow  water  approximation  used  here  as  a  basis  for 
the  theory  is  not  sufficiently  accurate  for  a  floating  beam 
whose  length  is  too  much  less  than  the  wave  length,  and  we  are 
now  trying  to  learn  something  about  the  magnitude  of  the  errors 
committed  in  using  the  shallow  water  theory  in  such  cases, 
Nevertheless,  it  does  seem  possible  to  design  floating  break- 
waters of  reasonable  length  which  would  be  effective  at  a  given 
wave  length.   Perhaps  it  is  not  too  far-fetched  to  imagine  that 


It  might  not  be  amiss  to  consider  the  physical  reason  why  it 
is  possible  that  a  beam  of  finite  length  could  be  more  effec- 
tive as  a  breakwater  than  a  beam  of  infinite  length.   Such  a 
phenomenon  comes  about,  of  course,  through  multiple  reflec- 
tions that  take  place  at  the  ends  of  the  beam.   Apparently 
the  phases  sometimes  arrange  themselves  in  the  course  of  these 
complicated  interactions  in  such  a  way  as  to  yield  a  small 
amplitude  for  the  transmitted  wave.   That  such  a  process  might 
well  be  sensitive  to  small  changes  in  the  parameters,  as  is 
noted  in  the  discussion,  cannot  be  wondered  at. 
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sections  could  be  added  to  or  taken  away  froni  the  breakwater  In 
accordance  with  changing  conditions. 

Another  consequence  of  the  theory--whlch  is  also  obvious  on 
general  grovinds--ls  that  there  Is  always  the  chance  of  creating 
a  large  standing  wave  between  the  shore  and  the  breakwater  because 
of  reflection  from  the  shore,  unless  the  waves  break  at  the  shore; 
this  effect  Is  perhaps  not  Important  If  the  main  Interest  Is  In 
breakwaters  off  beaches  of  not  too  large  slope  since  breaking  at 
the  shore  line  then  always  occurs.   (The  theory  developed  here 
could  be  extended  to  cases  In  which  the  shore  reflects  all  of 
the  incoming  energy,  it  might  be  noted.)   In  principle,  the 
calculation  of  the  deflection  curve  of  the  structure,  and  hence 
also  of  the  bending  stresses  in  it,  as  given  by  the  theory  is 
straightforward,  but  it  is  very  tedious;  consequently  only  the 
reflection  coefficients  have  been  calculated. 


1,   Theory  for  Floating  Obstacles  In  Shallow  Water. 

By  J.  J.  Stoker 

We  suppose  the  water  to  fill  a  region  lying  above  a  fixed 
horizontal  surface  (the  bottom)   y  =  -h,  and  beneath  a  surface 
y  =  Y(x,z;t),  the  motion  of  which  is  for  the  time  being  sup- 
posed known  (cf.  fig.  1.1).   The  y-axis  is  taken  vertically 


C 1 


y  =  -h 
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f-.. 
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>   Y(  X, s ; t ) 
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fig.  1.1 

upward,  and  the  x,z-plane  is  horizontal.   The  upper  surface  of 
the  water  given  by  y  =  Y(x,z;t)   will  consist  partly  of  the 
free  surface  (to  be  determined  by  tne  condition  that  the  pres- 
sure vanishes  there)  and  partly  of  the  surface  of  immersed 
bodies;  it  is,  however,  not  necessary  to  specify  more  about 
this  surface  for  the  present  except  that  it  should  represent 
always  a  small  displacement  from  a  rest  position  of  equilibrium 
of  the  combined  system  consisting  of  vmter  and  immersed  bodies. 

We  assume  that  the  motion  is  irrotatlonal  and  hence  that  a 
velocity  potential  £(x,y,z;t)  exists  such  that  the  velocity  q' 
of  the  water  particles  is  given  by   q  =  grad  ^  ;  this  function 
satisfies  the  Laplace  equation 


10 

(1.1)  ^   +  $■   +  ^   =0 
—XX   — yy   — zz 

In  the  space  filled  by  the  water.   The  letter  subscripts  here, 
and  in  what  follows,  denote  differentiations.   Since  we  deal 
with  small  motions  from  the  state  of  rest,  we  may  linearize  the 
basic  hydrodynamical  theory  and  ignore  all  quadratic  and  higher 
order  terms  in  ^  and  its  derivatives.   In  addition,  it  is  then 
legitimate  to  assume  that  all  boundary  conditions  at  the  upper 
surface  of  the  water  are  to  be  satisfied  at  the  equilibrium 
position;  this  position  is  supposed  given  by 

(1.2)  y  =  ^  (x,z) 

In  effect,  we  are  dealing  with  problems  which  concern  small 
oscillations  in  the  neighborhood  of  a  rest  position  of  equi- 
librium.  (The  bar  over  the  quantity  "f^  points  to  the  fact  that 
~  could  also  be  interpreted  as  the  average  position  of  the  water 
in  the  Important  special  case  in  which  the  motion  is  a  simple 
harmonic  motion  in  the  time.) 

At  points  on  the  free  surface,  rf  will  of  course  be  zero, 
and  this  in  turn  means  in  view  of  (1.2)  that  the  x,z-plane  is 
taken  in  the  undisturbed  position  of  the  free  surface.  Under 
the  immersed  bodies  the  value  of  ^^  will  be  fixed  by  the  static 
equilibrium  position  of  the  given  bodies.  Thus  ^  is  in  all 
cases  a  given  function  of  x  and  z;  for  a  floating  rigid  body, 
for  example,  it  would  be  determined  by  Archimedes'  law. 

The  condition  to  be  satisfied  at  the  upper  surface  is  the 
kinematic  condition: 

(1.3)  Ix  Y^  ^  3l^  Y^  -5y  +  Y^  =  0    , 

which  states  that  a  particle  once  on  the  surface  remains  on  it. 
At  the  bottom  surface,  of  course,  the  condition  to  be  satisfied 
is 


(1.3)-L  $y  =  ° 
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as  one  sees  on  physical  £;;rounds  as  well  as  from  (1.3)  applied 
to  the  case  Y  =  constant.   In  addition,  Bernoulli's  law  for 
determining  the  pressure  at  the  upper  surface  takes  the  form 


(1.1+)  E  +  1^  +  gy  =  0   for 


P 

when  quadratic  terms  In  the  velocity  components  are  neglected. 
We  now  write  the  equation  of  the  moving  upper  surface  In  the 
form 

(1.5)  y  =  Y(x,z;t)  =  -^(x,z)  -f-  7(x,z;t) 

and  assume  In  accordance  with  out  statement  above  that  'C(x,z;t) 
represents  a  small  vertical  displacement  from  the  equilibrium 
position  given  by  y  =  ^  .   Upon  insertion  In  (1.3)  and  (l.lj-) 
we  find  upon  ignoring  quadratic  terms  in  V   and  ^  and  their 
derivatives : 

(1.6)  ^rT  +^  "0     _^  +>0q^ 
'     -X  Lx   -z  iz   -y    ^t     / 

^   at   y  =  ■^(x,z) 

(1.7)  ^  -^  It  ""  S^T  ■"'^^  =  °  J 

as  boundary  conditions  to  be  satisfied  at  y  =  ^(x,z)   in 
accordance  with  remarks  made  above. 

At  points  corresponding  to  a  free  surface  where   p  =  0 
we  would  have,  for  example,  ^=0      and  hence 

(1.8)  -};y+<?t  =  o^ 

y  at   y  =  0 

(1.9)  2t  +  g'J  =  0  ,' 

These  are  the  usual  free  surface  conditions  (cf.  Lamb,  p.36ii). 
Another  special  case  might  be  that  in  which  the  motion  of  a 
portion  of  the  upper  surface  is  prescribed,  i.e.  <)(x,z;t)  as 
well  as  V   would  be  presumed  known;  in  such  a  case  the  condition 
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(1.6)  alone  would  suffice  as  a  boundary  condition  for  the  har- 
monic function  ^.   In  the  problems  to  be  treated  here,  however, 
we  do  not  wish  to  assume  that  the  motion  of  the  Immersed  body 
is  known  in  advance;  rather,  it  is  to  be  determined  by  the  inter- 
action with  the  water  which  exerts  a  pressure  p(x,z;t)  on  it  in 
accordance  with  (1.7).   Thus  the  exact  formulation  of  our  prob- 
lems would  require  the  determination  of  a  harmonic  function 
^(x,y,z;t)  in  the  space  between  y  =  -h  and   y  =  ^  which 
satisfies  the  conditions  (1.6)  and  (1.7)  at  the  upper  surface 
(in  particular  the  conditions  (1.8)  and  (1.9)  on  the  free 
surface)  and  (1.3),  at  the  bottom.   Additional  conditions  where 
immersed  bodies  occur  (that  must  be  obtained  from  the  appropriate 
dynamical  conditions  for  such  bodies)  would  be  necessary  to 
detennlne  the  pressure  p,  which  provides  the  "coupling"  between 
the  water  on  the  one  hand  and  the  immersed  bodies  on  the  other. 
Finally,  appropriate  initial  conditions  for  the  water  and  the 
immersed  bodies  at  the  initial  instant  would  be  needed  if  one 
were  to  study  non-steady  motions,  or — as  will  be  the  case  here-- 
conditions  at  oo  of  the  radiation  type  are  needed  if  simple 
harmonic  motions  (that  is,  steady  vibrations  instead  of  tran- 
sients) are  studied.   It  need  hardly  be  said  that  the  diffi- 
culties of  carrying  out  the  solutions  of  such  problems  are  very 
great  indeed — so  much  so  that  we  turn  to  an  approximate  theory 
which  is  based  on  the  assumption  that  the  depth  of  the  water  is 
sufficiently  small  and  that  the  immersed  bodies  are  rather  flat."" 

In  the  derivation  of  the  approximate  theory  we  start  from 
the  potential  equation  (1.1)  for  ^  and  Integrate  it  with  respect 
to  y  from  the  bottom  to  the  equilibrium  position"'*"  of  the  top 
surface   y  =  ^(x,z)   to  obtain: 


In  the  course  of  the  derivation  the  terms  neglected  are 
given  explicitly  so  that  a  precise  statement  about  them 
can  be  made. 

One  s_ees  readily  that  carryi_ng  out  the  integrations  to 
y  =  ri   rather  than  to  y  =  ^  +7   yields  the  same  results 
within  terms  of  second  order. 


','  I  f' 
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f- 
(1.10)     !  '  I       dy  =  -    (f   +  %^)    dy  =  | 

=  -  J_  i  '  ^  dy  -  A  (   (f)   dy  +  ^  (F  +  7   ^„ 
dx  (  ,  — X  -^   ©z  /  ,  — z  "^    IX  —X    <z  — z 
J  -n  1/  -n 

Here,  and  in  vhat  follovjs,  a  bar  over  the  quantity  ^  me  ens  that 
It  is  to  be  evaluated  at  the  equlllbriura  position  of  the  upper 
surface  of  the  water,  i.e.  for  y  =  ^{7i,z).      Prom  the  kinematic 
surface  condition  (1.6)  we  have  therefore  (duo  regard  being  paid 
to  the  fact  that  a  bar  should  nov;  be  put  over  ^  in  (1.6)): 

This  condition — really  a  continuity  condition--expresses  the 
fact  that  the  water  is  incompressible.   Consider  next  the  result 
of  integrating  by  parts  the  right  hand  side  of  (1.11);  in 
particular: 


3r  dy  =  <i^  +h^  -/ 


(1.12)     I  '  ar_  dy  =  TT  r  +  h  J  -  /  '  y  f _  dy 
v'-h  "^       '  ~^     -    J-h    ^ 


in  which  the  bar  under  a  symbol  means  that  its  value  Is  taken 
for  y  =  -h.   Since  we  have 


J 


If) 

'h^   dy  =  hS)'  -hcf 


(1.13) 

we  may  eliminate  ^  from  (1.12)  to  obtain: 

(l.lii)     !    f^  dy  =  (^  +  h)  ^^  -  /    (h  +  y)  (F__  dy 
w  -h  "^  ~^  k1  -h  ^y 

Indeed,  we  have  qalte  generally  for  any  function  P(x,y,z,t)  the 
result : 

( ^  -       -    r  ^ 

(I.IU),        F  dy  =  (7  +  h)  F  -      (h  +  y)  F  dy 

^    J-h  *  j-h  y 
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Makinei;  use  of  the  analog'ous  expression  for  the  Integral  of  ]| 
we  obtain  from  (1.11)  the  relation 

(1.15)     v?^  =  -r(Y  ^  h)  ?j^  -  r(v  ^  h)  1^]^  -  \^  J^ 

in  which 

(1.16)  I  =  j  ■  (h  +  y)  ?   dy,    J   =    j         (h  +  y)  ^   dy. 

j-h  ^y  j-h  ^y 

In  addition,  we  have  from  (1.10)  in  combination  v;ith  (1.14)  the 
condition: 

(1.17)  '^      =  -(^  +  h)^J        +  "3)  ]    ■*-  1      +  J  , 
— y     {      —XX   — zz     X    z' 

as  one  can  readily  verify. 

Up  to  this  point  we  have  made  no  aPDroximations  other  than 
those  arising  from  linearizing  for  small  oscillations.   The  es- 
sential step  in  obtaining  our  approximate  theory  is  now  taken 
in  neglecting  the  terms  I   and  J  .   This  in  turn  is  justified 
if  it  is  assumed  that  certain  second  and  third  derivatives  of  ^ 
are  bounded  when  h  is  small  and  that  v^  and  its  first  deriva- 
tives are  small  of  the  same  order  as  h:   one  sees  that  the 

terms    I      and   J      in  the   riiiiht  hand   sides    of    (1.15)    and    (1.17) 
X  z  2  ^^ 

are  then  of  order  h  v/hile  the  remaining  terms  are  of  order  h. 
Under  the  free  surface  in  the  case  of  a  simple  harmonic  oscil- 
lation one  can  show  that  this  approximation  requires  the  depth 
to  be  small  in  comparison  with  the  wave  length,  as  we  shall  see 
in  a  later  section. 

Upon  differentiating  the  relation  (1.7)  for  the  pressure 
at  the  upper  surface  of  the  water  v/i  th  respect  to  t  (again  not- 
ing that  a  bar  should  be  placed  over  the  term  ^.  in  (1.7))  and 
using  (1.15),  we  find  the  equation 

after  dropping  the  terms  I^  and  J  .   This  is  the  basic 

X.         z 
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differential  equation  for  the  function  ^(x,z)  which  holds 
everywhere  on  the  upner  surface  of   the  water.   In  particular, 
v;e  have  at  the  free  surface  where   p  =  0  and  v?  =  0  the 


equation 


(1.19) 


—XX   -zz   gh  -'-tt 


that  is,  the  linear  wave  equation  in  the  two  space  variables 
x,z  and  the  tiTie  t.   As  a  consequence,  all  disturbances  prop- 
agate with  the  constant  speed   c  =  Vgh  ,  as  is  well  known  for 
this  equation. 

If  there  is  an  immersed  object  in  the  v/ater,  the  equation 
(1.19)  holds  everywhere  in  the  x,z-plane  exterior  to  the  curve 
C  which  defines  the  water  line  on  the  immersed  body  in  its 
equilibrium  position.   The  curve  C  is  supposed  given  by  the 
equations 


(1.20) 


x=x(s),    z=z(s) 


in  terms  of  a  parameter  s.   We  must  have  boundary,  or  perhaps 
it  is  better  to  say,  transition  conditions  at  the  curve  C 
which  connect  the  solutions  of  (1.19)  in  the  exterior  of  C  in 
an  appropriate  manner  with  the  motion  of  the  water  under  the 
immersed  body.   Reasonable  conditions  for  this  purpose  can  be 
obtained  from  the  laws  of  conservation  of  mass  and  energy. 
Consider  an  element  of  length  ds  of  the  curve  C  representing 
the  water  line  of  the  immersed  body  (cf.  fig.  1.2).   The 


\ 


Z 


dsV 


/< 


C(x(s),y(s)) 


fig.  1.2 
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expression 

■:     Tv         ] 

^Pj_j,dy>ds 

represents  the  mass  flux  through  a  vertical  strip  having  normal 

n  and  extending  from  the  bottom  to  the  top  of  the  water.   From 

(1.1  1)t  applied  for  P  =  ^  we  have 
±  — n 

r^  _       -       r^ 

(1.21)   p     1^  dy  =  p(v/  +  h)  ^^  -  p  /    (h  +  y)  J        d.y 
J  -h  -'  -h 

-■  p(v^  +  h)  1^ 

2 
where  the  second  term  is  ignored  because  it  is  of  order  h  . 

Thus  it  v/ould  be  reasonable  to  require  that  (v^  +  h)  ^  should 
be  continuous  on  C  since  this  is  the  same  as  requiring  that  the 
mass  of  the  water  is  conserved  within  terms  of  the  order  re- 
tained otherwise  in  our  theory.   For  the  flux  of  energy  across 
a  vertical  strip  with  normal  n  we  have 

'        (     r^  r?  ) 

p  ^  dy  ^  ds  ^  -P  :   It-  I^  dy  -  gP    '  y  1^  dy  \   ds 
-h        ;      I       J-h     ^      ^  v/'-h    "    ) 

uDon  making  use  of  the  Bernoulli  law  (1.4)  for  the  pressure  n . 
Once  more  we  may  ignore  the  second  term  in  the  brackets  since 
it  is  of  order  h^.   Upon  apolying  (1.14)-j^  with  F  =  ^^  ^   and 
again  ignoring  a  term  of  order  h  we  find 

{  rv  ■        c    -        -  -  ) 

(l.S.-^)     )    \        P  1„  dyf  ds  =  :-p(>7  +  h)  f^   I4.  Ms . 

{  J-h        )      ''  n   r  , 

Since  we  have  alreadv  required  that  (v;  +  h)  ^   should  be 

'      — n      — 

continuous,  we  see  that  the  additional  requirement,  ^.  contin- 
uous, ensures  the  continuity  of  the  energy  flux. 

ks    transition  conditions  on  the  curve  C  delinitinp;  the 
Immersed  body  we  have  therefore 

(1.24)  (7  +  h)  ^  ,  ?,    continuous  on  C. 
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Of  course.  If  v^  is  continuous  (e.g.  if  the  sides  of  the  body  in 
contact  with  the  water  do  not  extend  vertically  below  the  un- 
disturbed free  surface)  it  follows  that  1  would  then  be  con- 

— n 

tinuous . 

In  order  to  make  further  progress  it  is  necessary  to  spe- 
cify the  properties  of  the  immersed  body.   However,  v;e  have 
succeeded  in  obtaining  the  equation  (1.18)  which  is  generally 
valid  and  of  basic  importance  for  our  theory  together  with  the 
transition  conditions  (1.24)  valid  at  the  edge  of  immersed 
bodies;  in  particular,  vie   have  the  definitive  equation  for  the 
free  surface  itself  in  the  form  of  the  linear  wave  equation 
(1.19).   The  idea  behind  this  method  of  approximation  is  to  get 
rid  of  the  depth  variable  by  an  integration  over  the  depth,  so 
that  the  problems  then  are  considered  only  in  the  x,z -plane. 
At  the  same  time,  the  problems  are  no  longer  problems  in  poten- 
tial theory  in  three  space  variables,  but  rather  problems  in- 
volving the  wave  equation  with  only  two  space  variables,  and 
hence  they  are  more  onen  to  attack  by  known  methods .   How  this 
comes  about  v/ill  be  seen  in  special  cases  in  the  next  tv;o  sec- 
tions . 


■  >\^\    r.ni   ^c. 


f.f 


1: 


*,  f!  ^i . 


•<  -'  '     .-> . 
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2 .   The  Effect  of  a  Floating;  Rigid  Body  on  SJTiple 
Harmonic  VJaves  .  " 

By  B.  Fleishman 

We  consider  only  two-dimensional  motions  of  a  rigid  body 
in  the  form  of  a  thin  plane  slab  (cf.  fig.  2.1).   The  ends  of 


^y 


x=-a 


x=a 


h 


fig.  2.1 

the  slab  are  at  x  =  +  a.   In  accordance  with  the  theory  pre- 
sented above  we  must  determine  the  siirface  value  ^(x,t)  and 
the  displacement  >y(x,t)  of  the  board  from  the  differential 
equations  (cf.  (1.15),  (1.19)): 


(2.1) 


—XX 


W^  %t' 


>   a 


(2.2) 


»?t  =  -h  1 


XX 


<  a 


We  have  dropped  the  bar  over  the  quantity  ^.  Tie   have  also  as- 
siomed  that  \^(x)  for  jxl  <  a,  the  rest  position  of  equilibriiam 
of  the  board,  is  zero;  this  is  an  approximation  that  is  justi- 
fied because  we  assume  that  the  board  is  so  light  that  it  does 
not  sink  appreciably  below  the  v/ater  surface  when  in  equili- 
brium.  (This  assumption  is  by  no  means  necessary--it  would  not 


"These  cases  have  been  treated  by  F.  John  in  the  paper  cited  in 
the  introduction. 
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be  difficult  to  deal  with  the  problem  if  this  simplifying  as- 
sumption were  not  made.) 

Since  "v?  is  zero,  it  follows  (cf.  (1.24))  that  the  transi- 
tion conditions  at  the  ends  of  the  board  are 

(S.r^)  ^  »  ^4-    continuous  at    x  =  +  a. 

We  are  interested  in  the  problem  of  the  effectiveness  of 
the  floating  board  as  a  barrier  to  a  train  of  waves  coming 
from  the  right  (x  =  +oo).   The  equation  (2.1)  has  as  its  gen- 
eral solution 

?(x,t)  =  F(x  -  ct)  +  G(x  +  ct),    c  =  ^. 

in  terms  of  two  arbitrary  functions  F  and  G  (as  one  can 
readily  verify)  which  clearly  represent  a  superposition  of  two 
Drogressing  waves  moving  to  the  right  and  to  the  left,  respec- 
tively, with  the  speed  \/gh  .   It  is  natural,  in  our  present 
problem,  to  expect  that  for  x  >  a   there  would  exist  in  gen- 
eral both  an  inconing  and  an  outgoing  wave  because  of  reflec- 
tion from  the  barrier,  while  for  x  <  a  we  would  prescribe 
only  a  wave  s:oing  outv>'ard  (i.e.  to  the  left).   vVe  shall  see 
that  these  qualitative  requirements  lead  to  a  unique  solution 
of  our  problem. 

We  consider  only  simple  harmonic  waves;  it  is  thus  natural 
to  write 

(2.4)  |(x,t)  =  i(x)  e^^^,       !x!  >  a, 

(2.5)  v^(x,t)  =  v(x)  e^^,    -a  <  X  <  a, 

with  the  stipulation  that  the  real  part  is  to  be  taken  at  the 
end.   (It  is  necessary  also  to  permit  {b(x)  and  v(x)  to  be 
complex-valued  functions  of  the  real  variable  x.)   The  condi- 
tions (2.1)  and  (2.2)  now  become 
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(2.6)  d_|  ^  a_  ^  ^  Q^    1^1  ^ 


dx^   -^h 


(2.7)  d_i  ^  i2  ^  ^  Q^    1x1  <  a. 
The  equation  (2.6)  has  as  general  solution 

(2.8)  (l>(x)  =  Ae^^  +  Be"^^^ 
with  k  given  by 

(2.9)  k  =  (T/Zgh    . 
For  ^(x,t)    we   have    therefore 

(2.10)  ^(x,t)    =  ^.e^(^-^^^    +  Be^('^^-'^^^\ 

the  first  term  representing  a  progressing  wave  moving  to  the 
right,  the  second  a  wave  moving  to  the  left.   In  our  problems 
we  prescribe  the  incoming  wave  from  the  right,  and  hence  for 
(b(x)  we  write 

(2.11)  hU)    =  Be^^  +  Re"^^^,    x  >  a, 

in  which  B  is  prescribed,  while  R--the  amplitude  of  the  re- 
flected wave  (more  precisely,  Ir'  is  its  amplitude ) --is  to  be 
determined.   At  the  left  we  write 

(2.12)  <l)(x)  =  Te^^'^, 

with  T--the  amplitude  of  the  transmitted  wave--to  be  deter- 
mined. 

To  complete  the  formulation  of  the  problem  it  is  necessary 
to  consider  the  dynamics  of  our  floating  rigid  body.   We  shall 
treat  two  cases:   a)  the  board  is  held  rigidly  fixed  in  a  hori- 
zontal position,  b)  the  board  floats  freely  in  the  water. 


1    t: 


.«.>Vf  Vv      bi'.VJT.  "■  ' 
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a )   Rigidly  Fixed  Board . 

If  the  board  Is  rigidly  fixed  we  have  v^(x,t)  =  0,  and 
hence  (cf.  (2.5))   v(x)  =  0.   It  follows  from  (2.7)  that  i 
vanishes  identically  under  the  board  and  hence  that  h{x)    is  a 
linear  function: 

(2.1."^)  <i)(x)  =  Yx  +  S,    -a  <  X  <  +a. 

Since  ^  (x,t)  furnishes  the  horizontal  velocity  component  of 
the  water,  it  follows  from  (2.1:^)  and  (2.4)  that  the  velocity 
under  the  board  is  given  by  ve    ,  i.e.  it  is  constant  every- 
where under  the  board  at  each  instant--a  not  unexpected  result. 

We  now  write  dov/n  the  transition  conditions  at  x  =  +a 
from  (2.3),  making  use  of  (2.4)  and  of  (2.11)  at   x  =  +a   and 
(2.12)  at   X  =  -a;  the  result  is: 

Be^^a  -.  Re-^^^  =  va  +  S 

Be^^^  -  Re-^^^  =  v/ik 

Te  ^^^  =  -ra  +  5 

Te-^^^  =  v/ik. 

Once  the  real  number  B--the  amplitude  of  the  Incoming  wave--has 
been  prescribed,  these  four  equations  serve  to  fix  the  con- 
stants R,  T,  •■',  and  6  and  hence  the  functions  ^(x,t)  and 
v9(x,t).   The  pressure  under  the  board  can  then  be  determined 
(cf.  the  expression  (1.4)  for  Bernoulli's  law)  from 


(2.15)  p(x,t)  =  -p  I  =  -pla  Mx). 


lat 


(Observe  that  the  quantity  Y  in  (1.4)  Is  zero  In  the  present 
case . ) 

In  terms  of  the  dlmensionless  parameter 

(2.16)  e  =  2a/A  , 
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the  ratio  of  the  lenp;th  of  the  board  to  the  wave  length  A  (cf 
(2.10))  on  the  free  surface  given  by 


(2.17) 


A  =  2Tr/^, 


the   solution  of    (2.14)    is 


P    _  errlBe 


2eTri 


errl    +   1 


Be 


29Tri 


)TTi     +    1 


(2.18) 


V    = 


em      Be^^^ 

a     *  etri    +   1 


8   =  Be 


errl 


The  reflection  and  transmission  coefficients  are  obtained 
at  once: 


(2.19) 


^r  = 


^t  = 


R 


err 


|B|    J  2"2 

Vl  +  e^rr"^ 


B 


Vl  +  e^TT 


2„2 


They  depend  only  upon  the  ratio  9  =  2a/A  ,  as  one  would  expect 

p     2 
They  also  satisfy  the  relation  C   +  C.  =  1,  as  they  should: 

this  is  an  expression  of  the  fact  that  the  inccning  and  outgo- 
ing energies  are  the  same*   The  following  table  gives  a  few 
specific  values  for  these  coefficients: 


! 

9        1        C.                 CI 
'          t                 r       1 

0.5 
1.0 
2.0 

.54 
.30 
.157 

.85 
.95 
.936 

I . 


< .. .. 


'1  >  ■ 


^■'-      -^^    fr^i  ^-f  f  <;'!i?in*t.^ 


!         i'  / 
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Thus  a  fixed  board  whose  length  is  half  the  incoming  wave 
length  has  the  effect  of  reducing  the  amplitude  behind  it  by 
about  50  percent  and  of  reflecting  about  72  percent  of  the 
incoming  energy. 

It  is  of  interest  to  study  the  pressure  variation  under 
the  board.   This  is  given  in  the  present  case  (cf.  (2.15))  by 

(2.20)  p(x,t)  =  -icrp(vx  +  5)e^'''^, 

the  real  part  only  to  be  taken.   Thus  the  pressure  varies 
linearly  in  x,  but  it  is  a  different  linear  fxinction  at  dif- 
ferent times  since  v  and  S  are  complex  constants  and  the  real 
part  only  is  to  be  taken.   The  result  of  taking  the  real  part 
of  the  right  hand  side  of  (2.20)  can  be  readily  put  in  the  form 

(2.21)  p(x,t)  =  P-,  (x)  cos  o-t  -  Pp(x)  sin  at 
with 

p,  (x)  =  apB(b-,  (x)  sin  r  +  bp(x)  cos  r) 

(2.22)  ,    r  =  9tt 
P2(x)  =  <TipB(bp(x)  sin  r  -  b^(x)  cos  r) 


and 


(2.2:^) 


1  +  r 


1  +  r 


We  have  assamed  in  making  these  calculations,  as  stated  above, 
that  B,  which  represents  the  amplitude  of  the  incoming  wave, 
is  a  real  number. 

In  fig.  2.2  the  results  of  computations  for  the  pressure 
distribution  for  time  intervals  of  l/4  cycle  over  the  full 
period  are  given  for  a  special  numerical  case  in  which  the 
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parameter  9  has  the  value   9=1,  i.e.  the  length  of  the  board 
Is  the  same  as  the  wave  length.   One  observes  that  the  pressure 
variation  is  greater  at  the  right  end  than  at  the  left,  which 
is  not  surprising  since  the  board  has  a  dampins;  effect  on  the 
waves.   One  observes  also  that  the  pressure  is  sometimes  less 
than  atmospheric. 

b)   Freely  Floating  Board. 

In  fig.  2.3  the  notation  for  the  present  case  is  indi- 
cated:  u(t),  v(t)  represent  the  coordinates  of  the  center  of 


^  X 


fig.  2.5 

gravity  of  the  board  in  the  displaced  position,  and  tO ( t )  the 
angular  displacement.   As  before,  we  consider  only  simple  har- 
monic oscillations  and  thus  take  u,  v,  and  Co   in  the  form 


(2.24) 


-   iCTt 
u  =  X  e    , 


-   icrt 
V  =  y  e 


CO  =  <>:>  e    , 


in  v/hich  x,  y,  and  oj  are  constants  representing  the  complex 
amplitudes  of  these  components  of  the  oscillation.   For  v^(x,t) 
we  have,  therefore 


(2.25) 


)^(x,t)  =  ^y  +  (x  -  x)co]  e 
=  (y  +  xZ3)  e^^^ 


iat 
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when  terms  of  first  order  in  x,  y,  and  60  only  are  considered. 
(The  horizontal  component  of  the  oscillation  is  thus  seen  to 
yield  only  a  second  order  effect.)   The  relation  (2.25)  now 
yields  (cf.  2.7): 

(2.26)  ^xx  "  ■  IT  ^y  ^  ^^^ 

in  which  h    is  the  complex  amplitude  of  the  velocity  potential 
5(x,t)  =  (|)(x)e  ^    .     Hence  (|)  is  the  following  cubic  polynomial: 

(2.27)  Mx)  =  -  ^  (^  -f  ^)  f  rx  +  8. 

Since  the  pressure  is  given  by  p  =  -p  ^^  -  Pg^  we  have  in  the 
present  case 

(2.28)  p(x)  =  [-icr  p(|)(x)  +  (y  +  xu))]  e^^*. 

The  transition  conditions  (2.?')  at   x  =  +a  now  lead,  in 
the  same  way  as  above  from  (2.11)  and  (2.12),  to  the  equations 

Ee^^^  +  Re-^^^  =  -  ^  (^  +  1|!)  ^  a-  ^  5 

—  2 
□  Ika   D  -ika     cr   /  co  a   ,  -„  \    1  ., 
Be     -  Re      =   '  m    ^— 2"  +  ^^  >  "  k  ^ 

(2.29)  X  _  ?   _  2 

^^-ika  ^  .  12  (.  i-L|l  ^1|_)  _  ar  +  5 

These  four  equations  are  not  sufficient  to  determine  the  six 
constants  R,  T,  co  ,  y,  y,  and  S.   We  must  make  use  of  the  dy- 
namics of  the  floating  rigid  body  for  this  purpose.   ''Ve  have 
the  equations  of  motion: 


(2.50)  F  =  Ma,    and   L  =  la 

at  our  disposal.   In  the  first  equation  P  and  M  are  the  total 


i     4     -'"1  -., 
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vertical  force  on  the  board,  and  its  mass;  and  a  is  the  vertical 
acceleration  of  its  center  of  gravity,  I  the  moment  of  inertia, 
and  oc  the  angular  acceleration.   These  dynamical  conditions 
then  yield  the  following  relations: 

(2.31)         I       p  dx  ^  M  V,    !   px  dx  =  165, 
J -a  J -a 

and  these  in  turn  lead  to  the  equations 


/- 


(2.32) 


(  ~     3      -   2  )  o_ 

Uio-pr-  ^(^^+^)+YX  +  S]  +y+xcl)  ^/dx   =   -Mcr^y, 

—     4     —   3  _   ^ 

-io-p[-  if(^^+^)    +  Yx^   +  Sx]    +  xy  +  x^o)  }  dx 

2  — 


In  the  first  equation  we  have  ignored  the  weight  of  the  board, 
since  it  is  balanced  by  the  hydrostatic  pressure.   The  equations 
(2.29)  and  (2.32)  now  determine  all  of  the  unknown  complex  am- 
plitudes . 

We  shall  give  the  results  for  the  case  of  a  purely  vertical 
motion  without  rotation  (though  the  general  case  presents  no 
special  difficulty);  i.e.  we  set  60  =  0,  and  solve  the  above 
equations.   The  result  is 


(2.33) 


y 

= 

BAD 
E 

S 

= 

E 

r 

= 

H 

T  =  BK^^Sm    (|)  .  II  (1), 


in  which 


.■    '.'",    '       <««■£;  ^t       /"I 


.n- 


-f    •  • 


.-^1 


A  =  e®^S  D  =  -  ^^^ 


E  =  -  f  G^A  .  2a  -  ^^  ^   219, 

a  p 


2  2       2  2 


F  =  -  g —  +  — :^ —  +  2a,    h  =  ■g;^  (err  -  1  j , 

a  p 


2a  \  =  ?J1 


A 


k 


From  the  expressions  for  R  and  T  In  equations  (2.19),  the 

reflection  and  transmission  coefficients,  C   and  C,  respec- 

I  R|  ItI 

tively,  can  be  calculated,  using  C   =  -r— -r  and  C,  = 


Again  one  finds  that   C^  +  C^  =  1. 

r    t 


r    B   "^^^   -^t    iBl 


Vi/e  wish  next  to  study  the  pressure  under  the  board.   For 
this  purpose  it  is  necessary  to  insert  the  function  h   given  by 
(2.2'7)  in  (2.28)  and  take  the  real  part.   The  calculation  is 
straightforward,  but  tedious,  and  gives  the  result: 

(2.54)       p(x,t)  =  p-,(x)  cos  at  -  Pg(x)  sin  at 

with 

(   Pt  (x)  =  g'^  P^^  [E  sin  r  +  F  cos  r] 


(2.35) 


p„(x)  =  ^^'^^'^  "g  [F  sin  r  -  E  cos  r] 
'^  C   +  D 


The  quantities  C,  D,  E,  F  are  given  by 


.-^2 


C  =  4  +  r^(Q  -  1) 


D  =  -  I  -  r(Q  -  3) 


E  =  ^2,-_^4(Q  _   ^)  _  ^2(Q  +  1)  _  2] 

(2.36)        +  5  ^r^(Q  -  1)^  +  4r^(Q  +  4)] 

+  [r^G.(Q  -  1)  +  r^Q(Q  +  1)  +20] 


=  5^[2r^  +  2r]  +  ^^r^(Q  -  1)^  +  4rQ  +  |] 
+  [-2r'^Q  -  2rO] 


with 


r  =  err 

0  -  1  -  -^ 

■^   —   T,         T 

^   a 

In  fig.  2.4  the  results  of  calculations  for  the  pressure 
distribution  in  a  numerical  case  are  given.   The  parameters 
were  chosen  as  follows: 

9=1,     h=lft,,     B=l,      a  =4  ft. 

M  =  0.3p/g,    a   =  4.46  rad./sec,    A  =  8  f  t . 

A  few  observations  might  be  made.   First  of  all,  we  note 
that  in  both  cases  the  pressure  variation  at  the  right  end 
(^  =  1),  where  the  incoming  wave  is  incident,  is  greater  than 
at  the  left  end.   This  is  to  be  expected,  since  the  barrier 
exercises  a  damping  effect  on  the  wave  going  under  it.   The 
pressure  distribution  in  the  case  of  the  floating  board  is 
quadratic  in  x,  in  contrast  with  the  case  of  the  fixed  board 
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in  which  the  distribution  of  pressiire  was  linear  in  x.   Next, 
we  note  that  the  pressure  variation  at  the  right  end  of  the 
stationary  board  is  greater  than  at  the  same  end  of  the  float- 
ing one;  this  too  might  be  expected  since  the  fixed  board  re- 
ceives the  full  impact  of  the  incident  wave,  while  the  floating 
one  yields  somewhat.   Finally,  we  see  that  at  the  left  end  the 
opposite  effect  occurs;   there  the  pressure  variation  under  the 
stationary  barrier  is  less  than  that  under  the  floating  one. 
This  is  not  surprising  either,  since  the  fixed  board  should 
damp  the  wave  more  successfully  than  the  movable  one. 
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3.   The  Effect  of  a  Floatin'^  Beam  vfith  Bending 
Flexibility  on  Simple  Har'nonic  Waves. 

By  L,  Weliczker 

The  beam  is  assiamed  to  extend  from  x  =  -^  to 
X  =  0   and,  as  in  the  above  section,  to  be  in  simple 
harmonic  motion  due  to  an  incoming  wave  from  x  =  +  oo  . 
The  basic  relations  for  ^(x,t)  on  the  free  surface,  and 
for  '^(x,t)  under  the  beam  are  the  same  as  before: 

(^•^^      ixx  =  iH^tt,  ^   >   0,    X   <   -l. 
(3.2)     71^  -  -h  ^^^,   -X<  X  <  0. 

We  assume  once  more  that  the  beam  sinks  very  little  below 
the  water  surface  when  in  oqa  ilibrium  (i.e.  very  little 
in  relation  to  the  depth  of  the  water) ,  so  that  the  coef- 
ficient of  T  in  (3«2)  can  be  taken  as  h  rather  than 
(h  +  y^)  (cf.  (1.15)),  and  also  the  transition  conditions 
at  the  ends  of  the  beam  are 


(3»3)     \       and  \        continuous  at   x  =  0,  x  =  -£,» 

After  writing 

(3.1^)     l{x,t)    =  (ti(x)e^'^S   )^(x,t)  -  v(x)e^^^ 

we    find,  as  before; 

_2 
(^•5)      *xx  ^  ih  't^  =  ^'   X  >  0,  X  <  -^ 

(^•^)      +XX  ^  TT  "^  =  °'  -^^   "^   <   ^' 

The  conditions  at   oo  have  the  effect  that  (cf.  (2.10) 
et  seq, ) : 

(3.7)      <t'(x)  =  Be'-^''  +  Re'^^"",   x  >  0, 
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(3.8)  i(x)  =  Te^^^,   X  <  -J,, 


with  k  =  crygh.   All  of  this  is  the  sane  as  for  the  previous 
cases.   We  turn  nov  to  the  conditions  which  result  from  the 
assumption  that  the  floating  body  is  a  beam. 

The  differential  equation  governing  small  transverse 
oscillations  of  a  beam  is 


(^•9)      EI  7^^^^^  ^  "^\tt  =  P' 


in  which  E   is  the  modulus  of  elasticity,   I   the  moment 
of  inertia  of  a  cross  section  of  unit  breadth  (or,  perhaps 
better,   EI   is  the  bending  stiffness  factor),   m   the  mass 
per  unit  area,  and   p   is  the  pressure,   '"'e  ignore  the 
vjeight  of  the  beam  and  at  the  same  time  disregard  the 
contribution  of  the  hydrostatic  pressure  term  in  p   corres- 
ponding to  the  equilibrium  position  of  the  beam i.e.   the 

pressure  here  is  that  due  entirely  to  the  dynamics  of  the 
situation.   Thus 


(3.10)     p  =  -  p  5^  -  ryr^r^ 


/         ■         I  \    ic^t 

=  (  -  lor  p(j)  -  pgv)e 

Insertion  of  this  relation  in  (3.9)  and  use  of  (3.6)  leads 
at  once  to  the  differential  equation  for   (}>(x): 

(3.11)    d!i  ^  pg  -  ma  d^ 

dx       EI    dx    Elh 


^  pg  -  ma^  dfi  ^  £g^  ^  Q^   -i<x<0, 


that  is  valid  under  the  beam.   The  case  of  greatest  importai  ce 

for  us that  of  a  floating  beam  used  as  a  breakwater leads 

obviously  to  the  boundary  conditions  for  the  beam  I'/hich 

correspond  to  free  ends,  i.e.  to  the  conditions  that  the 

shear  and  bending  moments  should  vanish  at  the  ends  of  the 

beam.   These  conditions  in  turn  mean  that  1)     and  T) 

ixx  (xxx 
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should  vanish  at  the  ends  of  the  beam,  and  from  (3*6)  we 
thus  have  for  cj)   the  boundary  conditions 

(3.12)     ^  =  ^=   0      at   X  =  0,  X  ^  -^. 
dx^   dx^ 

The  tr-^nsition  conditions  (3.3)  reqi  ire  that   (j)   and   (j) 
be  continuous  at   x  =  0,  x  =  -  ^,  Pnd  this,  in  view  of 
(3.7)  -nd  (3.B),  requires  that 

(<t)(0)  =  B  +  R   ,     (j)  (0)  =  ikB  -  ikR 


We  remark  once  more  that  the  constant   B   is  assumed  real, 
but  that   R   and  T  v.'ill  in  general  be  complex  constants, 
and  that  the  real  parts  of  J  and  7?     as  given  by  (3.1|.) 
are  to  be  taken  at  the  end. 

In  order  to  solve-  our  problem  vie  m.ust  solve  the  dif- 
ferential equation  (3.11)  subject  to  the  conditions  (3.12) 
and  (3»13)»   A  count  of  t  he  relations  available  to  determine 
the  solution  should  be  made:   The  general  solution  of  (3,11) 
contains  six  arbitrary  constants,  and  vje  wish  to  determine 
the  constants   R   and  T   (the  amplitudes  of  the  reflected 
and  transmitted  vjaves)  occurring  in  (3.13)  once  the  constant 
B   (the  amplitude  of  the  incoming  wave)  has  been  fixed.   In 
all  there  are  thus  eight  constants  to  be  found,  and  we  have 
in  (3,12)  and  (3.13)  eight  relations  to  determine  them.   Once 
these  constants  have  been  determint.d,  the  reflection  and 
transmission  coefficients  are  kno'-n,  and  the  deflection 
of  the  beam  can  be  found  from  (3.6).   The  maximum  bending 
stresses  in  the  beam  can  then  be  calculated  from  the  usual 
formula,   s  =  Mc/l,   with  M  =  EI  Ty,  .   and   c   the  distance 
from  the  neutral  axis  to  the  extreme  outer  fibres  of  the  beam. 
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In  principle,  therefore,  the  solution  of  the  problem  is 
straightforward.   Hovjever,  the  carrying  out  of  t  he  det-AL  Is 
in  the  case  of  the  bea.^  of  finite  length  is  very  tedious, 
involving  as  it  docs  a  system  of  eight  linear  equations 
for  eight  unknowns  v.'lth  complex  coefficients.   In  addition, 
one  must  determine  the  roots  of  a  sixth  degree   algebraic 
equation  in  order  t  o  find  the  general  solution  of  (3.11)« 
These  roots  are  in  general  complex  numbers  and  they  involve 
the  essential  parameters  of  the  mechanical  system.    Thus 
it  is  clear  that  a  discussion  of  the  behavior  of  the  system 
in  general  terms  viith  respect  to  arbitrary  values  of  the 
parameters  of  the  system  is  not  feasible,  and  one  must 
turn  rather  to  concrete  cases  in  vrhich  most  of  the  parameters 
have  been  given  specific  numerical  values.   The  results  of 
some  calculations  of  this  kind  have  already  been  given 
in  the  summary  at  the  beginning  of  this  paper,  though  only 
for  the  reflection  coefficient.   The  determination  of  the 
deflection  of  the  beam,  though  straightforward,  would  Involve 
a  great  deal  more  labor  than  is  recu  ired  to  determine  the 
reflection  coefficient. 

The  case  of  a  semi-infinite  beam i,e,  a  beam  extending 

from  X  =  0   to   x  =  -  co is  simpler  to  deal  I^ri  th  in  that 

the  conditions  in  the  second  line  of  (3.13)  fall  avray,  and 
the  conditions  (3.12)  at  x  =  -  .£  can  be  replaced  by  the 
requirement  that   (j)   be  bounded  at   x  =  -  oo.  The  number  of 
constants  to  be  fixed  then  reduces  to  four  Instead  of  eight, 
but  the  determination  of  the  deflection  of  the  beam  still 
remains  a  formidable  problem;  vje  shall  consider  this  case 
as  well  as  the  case  of  a  beam  of  finite  length. 

We  begin  the  program  indicated  with  a  discussion  of  the 
general  solution  of  the  differential  equation  (3.11).   Since 
it  is  a  linear  differential  eqi  atlon  with  constant  coefficients 
we  proceed  in  the  standard  fashion  by  setting   (|>  =  e'^  , 
inserting  in  (3.11),  to  find  for  K   the  eqa atlon 


vi  th 
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(3. 111.)    K^  +  n  K^  +  b  =  0 


(3.15) 


=  pp:-t^cr 


Elh 


This  is  a  cubic  equation  in  t   =  p,  which  happens  to  be  in 
the  standard  form  to  vrhich  the  Cardan  formule  for  the  roots 
of  a  cubic  appli-s  directly.   For  the  roots  0,   of  this  cubic 
equation  one  has  therefore 

p^  =  u  +  V 

2 
pp  =  eu  +  e  V 

2 
p,  =  e  u  +  ev 


(3.16) 


with  u   and   v   defined  by 

2 


(3.17)    u  = 


b 
2 


b 

T  ^  27 


3  i 
a^.2 


-  ^  +  (t-  +  ^) 


V 


b 
2 


3  t 


27) 


and   e   the  cube  root  of  unity: 
(3.18)    e  =  -^  ^3-V3  . 


The  constant   a   is  positive,  since   o' ,  the  f  reqa  ency  of  the 

incoming  wave,  is  so  small  in  the  cases  of  interest  in  practice 

2 
that  pg   is  much  larger  than  m&      ,      The  const.ant   b   is 

obviously  positive.   Consequently  the  root   p,   is  real  and 

negative  since   |u|  <  |v|   and  v   is  negative.   Thus  the 


roots  K". 


=  .  3^2, 


K-2  =  -Pi 


1/2 


are  pure  imaginary.   The 


quantities   pp   and   p^   are  complex  conjugates,  and  their 
square  roots  yield  two  pairs  of  complex  conjugates 


K3  = .  4^' 


n: 


k 


'  P 


1/2 


^5 


+  P 


1/2 

3  ^ 


^h  =  -  h 


1/2 
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For      Pp      and      i3        vje   have 


(3.19) 


=    -   i(u    +    V)    +    i^(u    -    V), 


(3.20)  P3  =   -  |(^  +  V)    -    i^(u   -   V). 


Thus   Pp   and   3^   both  have  positive  real  parts.   We  suppose 
the  roots  K'-,,  >f ,  ,   |r^,  ^,      to  be  numbered  so  that   f(, 
and   tfi;^   are  taken  to  have  positive  real  parts,  vrhile  K.  1 
and  l/[/      have  negative  real  parts.   The  general  solution 
of  (3.11)  thus  is 

K-iX      ^pX      i^^x      ^1  X 

(3.21)     <f'(x)  =  a,e    +  a^e  '   +  a-,e  -^   +  a,  e  ^ 

1^34 


K^x        k;^. 


X 


+  a^e  -^   +  a/C 


In  the    case   of  a    beam  covering    the  whole    surface   of  t  he 
vjater,    i.e.    extending   from      -   00    to      +   co  ,    the   condition   that 
d)      be   bounded   at      x   =   -   co    viould  require    that      a,   =  a,     =  a^- 

j>       k       s 

=  a,  =  0  since  the  exponentials  in  the  corresponding  terms 
have  non-vanishing  real  parts.   The  remaining  terms  yield 
progressing  vjaves  traveling  in  opposite  directions;  their 
wave  lengths  are  given  by  X  =  271/ K,  |  =  2T:/|h.',I   and  thus  by 

(3.22}     A  =  2VV  iu  +  v|  , 

with  u  +  V  defined  by  (3.I7).   The  vrave  length  and  frequency 
are  thus  connected  by  a  rather  complicated  relation,  and, 
unlike  the  case  of  vrave s  in  shallo'.'  vrater  vrith  no  immersed 
bodies  or  constraints  on  the  free  surface^  the  wave  length 
is  not  independent  of  the  frequency  and  the  wave  phenomena 
are  subject  to  dispersion. 
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In  the  case  of  a  bean  extending  from  the  origin  to 
-  00  v;hile  the  vjater  surface  is  free  for  x  >  0,   the 
boundedness  condition  for   (}>   at   -  oo  requires  that  -je  take 

)]/      have  negative  real  parts  and 


a/  =  0   since  ^fi   and  lu 


'Ir'^ 


K^x 


consequently  e  '^   and   e  ^   vjould  yield  exponentially 
unbounded  contributions  to  (|)   at   x  =  -  oo  .   "e  knovj  thrt   K-i 
and  >fp   are  pure  imaginary  v;ith  opposite  signs,  v;ith  IC , 
say,  negative  imaginary.   Since  no  progressing  v:ave  is  assumed 
to  come  from  the  left,  vie   must  then  take 


=  0 .   Thus  the 


term  a-,e 


K  X 


yields  the  transmitted  uave  and  the  terms 
Involving  a  and  a^  yield  disturbances  which  die  out 
exponentially  at  oo .  The  conditions  (3.12)  and  (3.13) 
at   X  =  0   novj  yield  the  following  four  linear  equation; 


(5.22) 


t^/a^  +  K^'a.  ^  If./'a^  =  0 


Jf-^^a^  +  rf^  a^  +  »%;a,  =  0 


^5  ■"    % 


=  B  +  R 


Jf^a,  +  dVa.,  +  K'r.a^  =  i^i^(3  -  R) 
for  the  constants   a-,  ,  a^,  a^  ,R .   For  the  amplitude  R   of 

X   o   *-' 

the  reflected  wave  one  finds 


(3.23) 


-1 

-1 

-1 

-1 

•^1 

«3 

^5 

ik 

/I 

'■1 

< 

'^a 

0 

R    - 

< 

o 

4 

0 

B 

-1 

-1 

-1 

-1 

^>i 

^^5 

ik 

^^1* 

K^ 

^^ 

0 

Kl 

Hf 

0 

\ 
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Even  in  this  relntively  simple  case  of  the  semi -Infinite 
bcain  the  reflection  coefficient  is  so  complicated  a  function 
of  the  parameters  (even  though  it  is  algebraic  in  then) 
that  it  seems  not  worthwhile  to  write  it  down  e::pllcitly. 
The  results  of  numerical  calculations  based  on  (3.23)  have 
already  been  given  in  Tables  A,  3,  C,  and  D  of  the  SLUx.iary 
at  the  beginning  of  this  paper. 

In  the  case  of  the  beam  of  finite  length  extending  from 
x  =  -X/to   x  =  0   the  eight  conditions  given  by  (3.12) 
and  (3.13)  must  be  satisfied  by  the  solution  (3.21)  of  the 
differential  equation  (3.11),  and  these  conditions  serve  to 
determine  the  six  constants  of  integration  and  the 
amplitudes   R   and  T   of  the  reflected  and  transmitted  v/aves. 
The  problem  thus  posed  is  qa  ite  straightforward  but  extremely 
tedious  as  it  involves  solving  eight  linear  equations  for 
eight  complex  constants.   The  amplitude   R   of  the  reflected 
v/avo,  for  example,  is  given  by  the  r^.tio  of  the  follov/ing 
eighth  order  determinants : 


K. 


-It? 


K- 


It; 


k: 


k: 


K' 


ir^e      K^e  K^e 


u ,  e      K^o  K^Q  fr^e     K  ,-e  Kr-e 


V 
-   1 

^1 


V 

-  1 


^^l 


fCgO 


-  1 


n. 


■^.J 


-  1 


lf= 


-M 


-^3^ 


B 


Zl 


0 


0 


0    0 

0    0 

0     0 
-10 


e  ^      0-0 


ik   0 

-iki 


-KqO 


0  -ikc"^'^' 
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in  which    y^    differs    frpn   the   dGtorminant    in   the   niunc-rator 
only   in   the    terra   in   the    seventh   column   and   fifth  rov;  v;hlch 
is      +1      rather    than      -    1    .      This    formula   was    used   to    calculate 
the   values    given  in   the    introduction   and   summary   at  the 
beginning  of  this    report. 
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